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Abstract 



We compute the helicity density matrix of p vector mesons produced in the 
two-body decays of polarized Xc2's in the framework of massless perturbative 
QCD. The Xc2's are either exclusively or inclusively produced in pp or pp 
interactions, via quark-antiquark annihilation or gluon fusion. Our results 
show unambiguous significant differences depending on the choice of the p 
distribution amplitudes and allow to discriminate between different proposed 
models. 

PACS numbers: 13.25. Gv, 12.38.Bx, 13.88.+e, 14.40. Gx 



1 



Typeset using REVTeX 



I. INTRODUCTION 



The general theoretical framework for the description of high energy exclusive processes 
within perturbative QCD has been widely discussed in the literature ^ and, although it has 
not yet reached the same level of reliability and variety of applications as the factorization 
theorem for the inclusive interactions, is expected to work at large momentum transfers and 
indeed several successful calculations have been performed so far in some simple cases . 
The numerical results depend strongly on the internal structure of the hadrons involved 
which, in the exclusive cases, are described by the hadron wave functions or distribution 
amplitudes, analogous to the distribution functions of the inclusive case. Perturbative QCD 
fixes the asymptotic form of these amplitudes ai Q'^ ^ oo and their general evolution, 
but in any realistic calculation they have to be taken as phenomenological quantities, to be 
experimentally determined via a set of physical information and then used in other processes. 
This procedure has led to the success of some distribution amplitudes motivated by QCD 
sum rules 

Some criticism has been advanced to the validity of the whole scheme, on the ground that 
it necessarily involves a contribution from soft non perturbative regions 0; such criticism 
has been answered by noticing that these soft contributions are suppressed by Sudakov form 
factors 10] . While this gives a better understanding and credibility to the theoretical model 
it weakens the agreement between theoretical computations, based on the QCD sum rule 
distribution amplitudes 0, and experiment; in fact these distribution amplitudes enhance 
the contribution of the soft end point regions, which are now known to be depressed by 
the Sudakov form factors. Moreover, some lattice estimates and further theoretical 
considerations |10,11| have suggested different shapes for the distribution amplitudes, closer 
to the asymptotic ones; the situation is then less clear and better phenomenological studies 
have to be performed in order to improve our knowledge of the properties of the hadrons 
and to reach the possibility of making genuine predictions. 

We consider here the Xc2 ^ P P decay process of polarized charmonium states created 
in pp or pp interactions and show how the observation of the polarization of the vector 
meson, via a measurement of its diagonal helicity density matrix elements, neatly depends 
on the p distribution amplitudes and helps in discriminating between different kinds of these 
quantities. 

The same process has already been studied, with different purposes, in two previous 



papers |]12| , p!3| , first in the exclusive channel pp ^ Xc2 P P |T2| and then in the inclusive 
one pp Xc2 + X pp + X [|13|. The goal there was that of assessing the significance 
of higher twist mass effects, rather than that of discussing the role of different distribution 
amplitudes: actually, it turned out that, independently of the chosen meson wave functions, 
only mass effects could give origin to non zero off-diagonal helicity density matrix elements, 
pi-i{p) or pifl{p). Any non zero measurement of even a single one of these non diagonal 
elements, signaling the violation of the helicity conservation rule, would be a clear signature 
of mass corrections. 

Here we compute, within massless perturbative QCD, the diagonal element pi,i(p) and 
study its dependence on the distribution amplitudes, ^{x)\ to be safe from mass corrections 
we exploit the results of Refs. and ||13[ and only consider kinematical regions in which 
such corrections are known to be small. The variations resulting from the choice of different 
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distribution amplitudes cannot then be ascribed to other reasons, 
we briefly recall the computations of Refs. 



12 and W 



In the next Section 
and show analytical expressions 
which are numerically evaluated in Section 3, which contains our results and where we also 
give some conclusive comments. We improve on the computations of Refs. |]12| and 
by explicitly introducing the evolution of the distribution amplitudes and the strong 
coupling constant as- 



II. COMPUTATION AND MEASUREMENT OF THE p HELICITY DENSITY 

MATRIX 

The full processes we are considering are either the exclusive, 

PP^Xc2^PP, (1) 

or the inclusive, 

PP^Xc2 + X pp + X , (2) 
production of a pair of p vector mesons with the subsequent decay 

p ^ TT TT , (3) 

and the quantity experimentally observed is the angular distribution of either one of the 
pions in the helicity rest frame of the decaying p . In the inclusive case (|]) the two initial 



particles need not be protons and might be other hadrons. 

The pion angular distribution depends on the spin state of the p via the elements of its 
helicity density matrix Pa,a'(p), 

3 

W{e, $) = — [po,o cos^ + (Pi,i - Pi,-i) sin^ 6 cos^ $ 
An 

+ + Pi,-i) sin^ e sin^ $ 

- V2{Re pi,o) sin 26 cos $] , (4) 

where B and $ are, respectively, the polar and azimuthal angles of the pion as it emerges 
from the decay of the p, in the p helicity rest frame. By integrating Eq. (H) over $ or one 
has respectively the polar and azimuthal distributions 

= ^[po,o + (Pi,i - Po,o) sin^ O] , (5) 

ly ($) = ^ [1 - 2pi,_i + 4pi,_i sin^ $] . (6) 
Zn 

Measurements of the above angular distributions yield direct information on Pa.a'(p)- 

On the other hand, the vector meson helicity density matrix p(p) can be computed in 
perturbative QCD 

P\,\'{p\0,(p) = ^ Ax,f,;MAy^^.M' PM,M'iXc2) (7) 
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where is the normahzation factor such that Tr[ p ] = 1, 



Ar(^,0)= ^A,m;m(^,0)AI,^;M'(^,0)Pm,M'. (8) 

X,^i,M,M' 

The Ax^^-m{(^, 0)'s are the hehcity amphtudes for the decay of a Xc2 at rest, with spin third 
component Jz = M, into two p vector mesons at angles 6, (p and n — 6,tt + (p and hehcities 
A and /i respectively, Xc2{Jz = M) p(A) p(p); Pm,m'{Xc2) is the spin density matrix of the 
decaying charmonium state. 



Explicit expressions of the helicity decay amplitudes Ax^^-m can be found in Ref. [|12 
computed within perturbative QCD in the general case in which (constituent) quark masses 
are taken into account; the massless case, which we consider here as explained in the Intro- 
duction, is simply obtained by setting e = in the expressions of Ref. ||T^. These decay 
amplitudes contain the p distribution amplitudes, (p{x). 

In order to compute the p density matrix via Eqs. (|^) and we still have to specify 
the spin state of the Xc2] to do so we treat separately the exclusive (P and inclusive 
production processes. The charmonium state is produced in the exclusive pp channel via 
quark- ant iquark annihilations and the helicity conservation rule of massless perturbative 
QCD entails Q: 

p^,^, = M^M' 

Pofo = Pt2 = P-2,^2 = (9) 
P^i = P-1,-1 = 1/2 . 

In the inclusive channel (|^) the Xc2 is produced dominantly via two gluon fusion and, 
assuming a non relativistic wave function for the charmonium state, one finds [l^ : 



Pl^ = pTi = = (10) 

pZ = P-2,-2 = 1/2 . 



By exploiting Eqs. ([7|)-([10|) and the analytical results of Ref. one derives for the 



helicity density matrix of the p meson, produced either via the exclusive process (|I|) or the 
inclusive one (0): 

Pa,a' = A ^ A' 

Po,o = 1 - 2pi,i (11) 
1 1 



Pi,i = P-1,-1 = 7T 



2 Uo.o 



2 



k1i,-i 



where the reduced amplitudes Ax y are given in Eqs. (2.10, 11) and (2.14, 15) of Ref. |T2 
and do not depend on the p production angles 9 and 0, but do depend on the distribution 
amplitudes. Eq. (|1^) has the same form both for exclusive and inclusive production of p, 
but the dependence on the production angle 9 is different in the two cases 
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1 + cos^ 

F'^"(0) = ^ (-lo) 

^ ' 1 + 6 cos2 e + cos4 e' ^ ^ 

Eq. (^) shows that in massless perturbative QCD only the diagonal elements of p{p) 
survive and their values, which depend on the p distribution amplitudes, will be discussed 
in the next Section; Eq. ([Tl|) also implies that the azimuthal angular distribution (|^) of 
the pions emitted by the p vector meson is fiat, iy($) = l/{2n), whereas the polar one @ 
depends on pi^i (or po,o) and can be rewritten as 

W{e) = ^ [pi,i + (1 - 3pi,i) cos^ 0] (14) 
= ^[2po,o + (l-3po,o)sin2e]. (15) 



III. NUMERICAL EVALUATION OF pi,i AND DEPENDENCE ON THE 

DISTRIBUTION AMPLITUDES 



We report here for convenience from Ref. |T2| the explicit expression of |Ao,o|/|^i,-i| 
which, together with the F{6) functions given in Eqs. (|T^) and (|l^), enters the computation 
of the helicity density matrix pi i according to Eq. (|TT|): 



I A 



0,0 




(16) 



where 



1 /■! 



^ as[xyMl]a,[il-x){l-y)M^^] 
xy{l - x)(l - y){2xy - x - y) 



'0,0 



1 /•! 

"32/0 ^^^y^^^^'^'^^^'^^^y^^l^ 

as[xyMl] as[{l-x){l-y)Ml] 
xy{l - x){l - y){2xy - x - y) 
{x - yf 



1 + 



2xy — X — y 



(18) 



In the above equations is the Xc2 mass, = min(z, 1 — z)Q |jT5| and (pi, are the 
p distribution amplitudes which, following Refs. |]1]J^, are in general assumed to be different 
for longitudinally (L, A = 0) or transversely (T, A = ±1) polarized vector mesons; fi and fx 
are the corresponding decay constants. 
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To test the dependence of the numerical values of pi^i on the distribution amplitudes 
we have performed the computations of the integrals ([T7|) and ([T8|) choosing two typical 
different sets of and ipx, i.e. the symmetric distribution amplitudes 



and the QCD sum rule ones 



ifiix) = '~Pt{x) = 6x(l - x) 



= 6x(l — x) 

X + (2a;-l) 



(19) 



'asiQlY 


2/3 1 







(20) 



V9T(a;, Ql) = 6x(l - x) 



-| 4/25 



- -Cf (2x-l) 





52/75 1 







(21) 



where /i| = 0.5 (GeV/c)^, /i^ = 0.25 (GeV/c)^ [0, and C{z) denotes Gegenbauer polyno- 
mials. In both cases we have fi = fr- 

The latter set includes explicitly the QCD Q"^ evolution of the distribution amplitudes 
|1.16] and the strong coupling constant has a smooth limited small behaviour [17 



127T 



25 In 



(Q2 + 4m2)/A2 



(22) 



with rrig = 0.5 GeV/c and A = 0.2 GeV/c. The effective gluon mass rrig, which avoids 
enhancements in the soft end point regions, models in a convenient way the effects of Sudakov 
form factors 0. 

Notice that Eqs. (^) and (|2T|), although they imply the usual asymptotic (Q^ — > oo) 
behaviour [llll 



asm 

asif^r) 



4/25 



(23) 



lead, at the moderate values involved in the processes considered here, to longitudinal 
and transverse distribution amplitudes which strongly differ in their x dependences; this 
feature strongly affects the numerical results. 

The set (|^), on the contrary, has the same x dependence both for (pi and (pr and no 
evolution; the same x dependence is the dominant feature and we have checked that 
by adding to ifx a dependent factor in agreement with Eq. (|23|) the numerical results 
for the helicity density matrix do not change significantly. This is the reason why we have 
only considered the two choices of (pi and ipT given above: they are representative of whole 
classes of distribution amplitudes. 

In Figs. 1 and 2 we show the resulting values of pi,i(p) as a function of the production 
angle 9 for the two choices of the distribution amplitudes (|T^) and (^ ^T|) and, respectively. 



6 



for the exclusive and inclusive production processes. In both cases the two sets of distribution 
amplitudes lead to clearly different results. Notice that pi,i(6') is symmetric around 9 = 7r/2, 
see Eqs. (0) and (|T3|). 

The values of pi^i fix the shape of the polar angle distribution (0) of the pion resulting 
from the p decay, which is shown in Figs. 3 and 4, respectively for p vector particles 
produced at = 50° in the exclusive case and 6 = 90° in the inclusive one. We have chosen 
such angles because in these kinematical regions the mass corrections are negligible and 
cannot affect our study of the distribution amplitudes; this can be seen from the results 
of Refs. ]T^ , p!B| (which, however, were derived with a fixed coupling constant and no 
evolution of the distribution amplitudes) and has also been explicitly verified with the present 
improvement of the calculations. Again, the two choices of distribution amplitudes, both in 
the exclusive and the inclusive channels, lead to very different pion angular distributions, 
whose observation would allow to discriminate between them. 

Our study shows that different distribution amphtudes give different results for the di- 
agonal elements of the p vector meson helicity density matrix; the p's are decay product of 
charmonium states Xc2, either exclusively or inclusively produced. A measurement of pi^i is 
difficult and might not allow a detailed comparison between similar distribution amplitudes, 
but certainly allows to discriminate symmetric distribution amplitudes, Eq. (|T^, or other 
distribution amplitudes with ip^ ~ ^pr, from those similar to the ones inspired by QCD 
sum rules, Eqs. (^, [2l| ). Even for such a qualitative information is worth attempting the 
measurements we propose here; the more we know about subtle hadronic properties and the 
more we gain in predictive power. 
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FIGURES 



FIG. 1. Values of pi^i{p) as a function of the p meson production angle 9, in the exclusive 
process (|l|). Solid curve: symmetric distribution amplitudes, Eq. (19); dashed curve: QCD sum 
rule distribution amplitudes, Eqs. ( pO| , [2l| ). 

FIG. 2. Values of pi^i{p) as a function of the p meson production angle 0, in the inclusive 
process Same notations for the curves as in Fig. 1. 

FIG. 3. Plot of the angular distribution, W{Q), of the vr emitted in the p decay. The p has 
been produced in the exclusive channel at an angle 9 = 50°. Same notations for the curves as in 
Fig. 1. 

FIG. 4. Plot of the angular distribution, W{Q), of the vr emitted in the p decay. The p has 
been produced in the inclusive channel at an angle 9 = 90°. Same notations for the curves as in 
Fig. 1. 
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